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Classical Gluon Production Amplitude for Nucleus—Nucleus Collisions: 
First Saturation Correction in the Projectile 
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We calculate the classical single-gluon production amplitude in nucleus-nucleus collisions in¬ 
cluding the first saturation correction in one of the nuclei (the projectile) while keeping multiple¬ 
rescattering (saturation) corrections to all orders in the other nucleus (the target). In our approxi¬ 
mation only two nucleons interact in the projectile nucleus: the single-gluon production amplitude 
we calculate is order-y® and is leading-order in the atomic number of the projectile, while resum¬ 
ming all order-one saturation corrections in the target nucleus. Our result is the first step towards 
obtaining an analytic expression for the first projectile saturation correction to the gluon production 
cross section in nucleus-nucleus collisions. 

PACS numbers: 12.38.-t, 12.38.Bx, 12.38.Cy 


I. INTRODUCTION 


Complete understanding of heavy ion collisions is impossible without a clear picture of the initial stages of the 
collision, i.e., the processes leading to the creation of quarks and gluons which later on thermalize forming quark- 
gluon plasma (QGP). The distribution of quarks and gluons produced in these early-time processes, known as the 
initial condition for the QGP formation, is the fundamental building block of heavy ion theory. Questions concerning 
thermalization of the quark-gluon medium and the determination of the initial conditions for the hydrodynamic 
evolution of the QGP can not be answered in a fully satisfactory manner without the qualitative and quantitative 
knowledge of the production mechanism for the initial-state quarks and gluons in a nuclear collision. 

In the framework of saturation physics (see the reviews and the book Q) the leading-order contribution to 
gluon production is given by the classical gluon fields (9l-ll3l| of the McLerran-Venugopalan (MV) model (l^ - [l^ . 
Glassical gluon fields in heavy ion collisions resum powers of of and aj. aT [13 , where Og is the strong coupling 
constant, while A\ and A2 are the atomic numbers of the two nuclei (henceforth referred to as the projectile and the 
target). Since the saturation scales squared of the two nuclei are proportional to these parameters, ^ and 

Qs2 ~ write down the quasi-classical single-gluon production cross section as 

da I ( QlAB±-b^) Ql2{b±) \ 

d^kd^B d^b Og y y 


where Bi_ is the impact parameter between the two nuclei, b± is the transverse position of the produced gluon with 
respect to the center of the target nucleus, while k± is the transverse momentum of the produced gluon with kp = |fcj_|. 
The expansion in the powers of o^aJ^^ and alA^^ corresponds to expansion in the powers of Qh/k'^ and Q^2/^t^ 



Note that due to projectile-target symmetry (resulting in Qgi -o- Qs2 symmetry) we have Cn,m = Cm,n- (The kp- 
dependence on the right-hand side of Eq. (ED also enters through powers of ln(fcT/A) where A is the infrared (IR) 
cutoff of each nucleon; the powers of In(fc7’/A) are not shown explicitly above. The coefficients Cn,m from Eq. ([ 2 |) are, 
in fact, polynomials in In(fc7’/A).) 

At the moment we do not have an analytic exp ression for the function f{Qli/k^,Q'^2/^T) (ED- This func¬ 

tion was extensively studied numerically in |I8l - l24| . Still it appears desirable to attain a better handle on the 
analytic form of this function. Apart from the general advantage of having an analytic solution, knowing the func¬ 
tion /(Qsi/fcxi Qs2/^t) should greatly facilitate the inclusion of small-a; evolution corrections [ 25 l - [ 33 | along with 
the running-coupling corrections |35l - l39l| into the gluon production cross section. These corrections are essential for 
realistic phenomenological applications. 
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Let us briefly summarize what is known about the function /(Q^i/^T’Qs 2 /^t)- The leading-order result (the 
coefficient Ci.i in Eq. @) was obtained in ilia, El, reproducing the earlier results of [H IH . The case of proton- 
nucleus [p + A) collisions, defined as the leading-order in term in the expansion on the right-hand side of Eq. ([2|), 
was solved in (see also [H, El), yielding the coefficients ci,„ (and, due to target-projectile symmetry, c„q) for 


any positive integer n. This is all we presently know analytically about the coefficients Cn,m- 

An ansatz for the full solution of the classical gluon production problem was proposed by one of the authors 


While consistent with our knowledge of 


in 0. A variational approach to the problem was attempted in 
coefficients ci,„, neither of these results can be verified further due to our lack of knowledge of the coefficients Cn,m 
for n,m >2. In phenomenological applications one often employs the A:T-factorization formula involving unintegrated 
gluon distributions : while the gluon production cross section in p -|- A collisions (the lowest-order in Qh terms 

in Eq. ([2])) does lead to the A:'r-factorization formula it is not clear whether fcT-factorization holds beyond 

the p + A approximation, again due to our lack of knowledge of Cn,m for n,m>2. Moreover, numerical simulations 


of the classical gluon production |18l - l24j for nucleus-nucleus (A -|- A) collisions appear to rule out the /cr-factorization 
ansatz, suggesting that this factorized result is only valid for p -I- A. 

The goal of the current project is to determine the coefficients C 2 ,n (for n > 2). As one can see from Eq. ([21), to 
obtain C 2 ,n coefficients one has to expand the gluon production cross section to the second order in Qsi keeping all 
orders of Q^ 2 - This means that one has to allow two nucleons in the projectile nucleus to participate in the interaction, 
while allowing all nucleons in the target nucleus to interact. Formally one can think of this as working in the regime 

where A^^ ~ 1 while A^^ <C 1: the goal is to calculate the O (o;^ = O A 


2/3 


correction to 


gluon production. Note that one still has Ai ^ 1 such that the 0{al) contribution to the cross section, where only 
one of the projectile nucleons interacts, is small due to it being suppressed by a power of A^ . 

As we will detail below, having two nucleons interact in the projectile nucleus results in the gluon production cross 
section consisting of two contributions (see the top two diagrams in Fig.|2]below): (i) each of the two nucleons interacts 
both in the amplitude and in the complex conjugate amplitude (or one nucleon interacts only in the amplitude and 
another nucleon interacts only in the complex conjugate amplitude), or (ii) one nucleon interacts both in the amplitude 
and in the complex conjugate amplitude while the other nucleon interacts only in the (complex conjugate) amplitude. 
In the case (i) the scattering amplitude is such that the contribution to the cross section is ~ \g^\^ ~ 

In the case (ii) the amplitude is O(g^), while the complex conjugate amplitude is 0(g), such that the cross section 
contribution is ~ ^ of. (As we will argue below, the case where the amplitude is 0(g'^) and the complex 

conjugate amplitude is O(g^) is included by using retarded Green functions in the O(g^) and O(g^) amplitudes in 
this quasi-classical calculation.) 

In the present paper we calculate the scattering amplitude in case (i). While its square would give a contribution 
to the desired gluon production cross section at order-Qfj^, the complete expression for the cross section can only 
be obtained if one includes the contribution of case (ii) as well. Calculation of the contribution (ii) would involve 
the O(g^) amplitude, which appears to be significantly more involved and is left for future work. Our calculation is 
performed in the light-cone gauge of the projectile, with the final results given in Eqs. (I36l) and (l40l) . 

The paper is structured as follows. In Sec. [IT] we set up the problem by first reproducing the lowest-order (p + A) 
gluon production calculation and then by outlining the main ingredients needed to complete the calculation of the 
coefficients C 2 ,n in Eq. The elements of the calculation, along with the main results, are presented in Sec. IIIII We 
conclude in Sec. lYl 


II. THE SETUP 


Consider the high-energy scattering of a projectile nucleus on a target nucleus. Defining the light-cone variables 
as = (v° ± v^)/^/2 with the a:^-axis being the collision axis, we choose the projectile to be moving along the x'^ 
light-cone direction, and the target moving along the x~ direction. All of our calculation will be performed in A’*' = 0 
light-cone gauge. Transverse plane vectors are denoted as El = (u^, u^), such that the full 4-vector is = (z;+, Ej_) 
in the (-1-,—,T) notation. 


A. Classical gluon production in the p -\- A approximation 


To introduce our formalism, let us begin by outlining the calculation of the single gluon production cross section 
in p -|- A collisions in the quasi-classical approximation. Working in the approximation where A^^^ ~ 1 and 
<C 1, in this Section we calculate the O a^A^^^/os = O Og A}^^ contribution. Since only one power of 
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A-^ is required we only need to include diagrams where one nucleon from the projectile interacts with the target. 
The diagrams contributing to the gluon production cross section in p + A are shown in Fig. [TJ The target nucleus 
moving along the x“-axis generates the shock wave, shown in this paper as a vertical band. Diagrammatically this 
vertical band represents multiple interactions with the field of the shock wave, which happen over a very short time 
interval around x~^ = 0. Since only one nucleon in the projectile is involved in the interaction, we model it with a 
single quark, shown in Fig. [U by a horizontal solid straight line. The spectator quarks are not shown explicitly in 
Fig. [T] (along with the spectator nucleons if the projectile is a nucleus): below, for other processes, we also do not 
show the spectators explicitly. The produced gluon can be emitted either before the interaction with the shock wave 
(left diagram in Fig.[T|) or after the interaction (right diagram in Fig.[T]). Emission during the passage of the projectile 
through the shock wave is suppressed by a power of center-of-mass energy. Gluon emission from within the shock 
wave is also suppressed in the A'^ = 0 light-cone gauge we are using. 






2± 


FIG. 1. Diagrams contributing to the single-gluon production amplitude in p + A collisions. 


Normalizing the incoming quark in the projectile as a Wilson line, we write the following expression for the scattering 
amplitude resulting from the diagrams in Fig. [T] 


A{z±,b±) = 2g 


<Pk 


■ ki 


kl 


jjab _ jjab 
•^-L 0_L 




(3) 


The amplitude in Eq. ([S]) includes a Fourier-transform from transverse momentum space to transverse coordinate 
space. It also involves the adjoint 


c|f= PexpL/..+rM-(A.-=o,=p 


ab 


(4) 


and fundamental 


Vr =Pexp<i (7 / A’^ {x~^,b =0,b±) 


(5) 


Wilson lines describing respectively the propagation of a high-energy gluon and quark through the target gluon field 
Here T“ and are correspondingly the adjoint and fundamental SU(fVc) generators. In arriving at Eq. Q we 
used the Fierz identity U°’^ = t°‘V when evaluating the right diagram in Fig.[T]in order to put the contributions 

of both diagrams into a similar form. 

Performing the Fourier transform in Eq. © we obtain 


A{z±,b±) = 


ig ■ {zx - b±) 


\z±-b±\^ 

The gluon production cross section is given by (see e.g. [j 
da 1 


Ug - U. 

•Z_L t) 


ab 




d^kp dy 2 (27r)^ 


J d'^zd^z' ^H(fj_,&j_)H*(f[,6j_)J>, 


( 6 ) 

(7) 


where the summation over colors and polarizations of the final-state particles along with the averaging over the 
quantum numbers of the initial-state particles are implicitly implied. The angle brackets (...) denote averaging in 
the target nucleus wave function. 
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Substituting the amplitude Eq. ([6]) into Eq. d?]) yields [41| 
da agCp 


d?kTdy 


d^z dPz' d'^b 


z± - bj 


z\ - bx 


\z±-b±\^ 1^1 

Sg{za., z'x_) — Scibi., -?!) — S'g(-?J-, &_l) + 1 


where the gluon dipole 5'-matrix is defined by 


Sail^A) = -pE_ (t,g cty) = ^ (ll' [U,, Ul] ) . 




In the quasi-classical MV/Glauber-Mueller (GM) approximation it is given by [M| 


SG{xx_,yi_) = exp 


- yifQlG 


x± + yj 


In 


x± - y±\ A_ 


( 8 ) 

(9) 


( 10 ) 


where (5 ^q(6_l) = 47ra^ T{b±) is the square of the gluon saturation scale with T{b±) the nuclear profile function and 
A the IR cutoff of each individual nucleon (A ^ Aqcd)- Eq. (fTOl) resums all the multiple rescatterings in the target 
nucleus: thus, in the MV model, it resums all-order saturation corrections in the target. 

Note a particular convenience of the form of the amplitude given in Eq. ©, with the fundamental Wilson line 
V placed to the left of the fundamental SU(A’c) generator in both terms there: when we square the amplitude, the 
fundamental Wilson line V is multiplied by its hermitean conjugate giving an identity. Below when calculating 
diagrams we will always cast them in the same form, which would make all fundamental Wilson lines vanish when we 
square the amplitude. 


B. First Saturation Correction in the Projectile: Diagram Types and Calculational Simplifications 

The goal of this project is to calculate the first projectile saturation correction to the cross section in Eq. (jS]). That 
is, we need to find the order-a^ correction to dH): this is the leading in Ai part of the order-a^ correction. To get 
the leading-Ai contribution, the correction must include an interaction with another nucleon in the projectile nucleus. 
Hence we see that we need to find order-a^ correction to dHI involving one other nucleon in the projectile. The main 
types of the diagrams we need to calculate are shown in Fig. [2j There by two horizontal solid straight lines we show 
two quarks from two different nucleons in the projectile nucleus, again suppressing the spectators. The diagrams 
in Fig. [2] represent the amplitude squared which contributes to the cross section: the solid vertical line denotes the 
final-state cut. The cross labels the measured gluon. 

The three types of projectile saturation corrections to gluon production shown in Fig. [2] are (i) the square of the 
order-g^ amplitude; (ii) the interference between the order-g® amplitude and the leading-order (oidei-g) amplitude 
from Eq. dSD; and (iii) the interference between the order-f;"^ amplitude and the oidei-g^ amplitude. (In this power 
counting we are assuming that the interaction with the shock wave is order-1, that is ^ 1, and are counting 

only powers of g arising from the vertices shown explicitly in Fig. [^) The order-a^ diagrams in which one of the 
nucleons is a spectator and does not participate in the interactions are suppressed by Aj^' and are neglected in 
our analysis. Note that diagram (iii) has two gluons in the final state: since we are calculating the single-inclusive 
production cross section, having more than the measured gluon in the final state is allowed. 

Below we calculate the order-amplitude that enters diagram (i) in Fig. [21 thus constructing the order-Og correction 
to the amplitude in Eq. (|6]) which is enhanced by the leading power of A\^^ . The order-g® amplitude from diagram 
(ii) is left for future work. 

The situation with the diagram (iii) is more subtle. First let us note that, since we are working in the classical MV 
model, the diagrams we consider for gluon production also correspond to the diag rams contributing to the classical 
gluon field with the two colliding nuclei providing the source current I, EE Ildl-IE. Hence one can think of these 
diagrams as graphically representing the solution of the Yang-Mills equations E, El • In such case each Feynman 
propagator would be replaced by the retarded Green function. While such a replacement is relatively straightforward 
for gluon propagators, it is less clear what this prescription means for quark propagators in the time-ordered picture 
we employ here: for instance, the solution of Yang-Mills equations can only have the produced gluon emitted by the 
projectile quark after the interaction with the shock wave (see e.g. El): whereas a naive application of perturbation 
theory allows one to draw a number of a priori non-zero graphs with the produced gluon emitted before the shock 
wave interaction (see Fig. |7| below). 
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(i) 




(Hi) 


FIG. 2. The diagrams representing different types of saturation corrections in the projectile nucleus to the gluon production 
cross section. 


To clarify this issue we need to perform a detailed diagrammatic analysis. First let us review how the gluon Feynman 
propagators become retarded Green functions. As an example consider the top two diagrams shown in Fig. [3] where 
regular Feynman propagators are implied for uncut lines. The difference between the top two diagrams in Fig. [3] is 
due to one of the vertices involving the gluon carrying momentum I being moved across the cut in the right diagram. 
The rest of the diagram is the same in both cases. Concentrating on the propagator of the gluon carrying momentum 
I and suppressing the rest of the diagrams’ contribution we see that adding the two graphs gives 




P+i€l+’ 


( 11 ) 



FIG. 3. An example of diagrams which add up to convert a gluon Feynman propagator into a retarded Green function. The 
bold arrow on the gluon line in the last diagram denotes a retarded gluon Green function, with the arrow’s direction pointing 
in the direction of light-cone time flow, that is, indicating that > x^ . 






































where 


6 


^M^(0 = 9tiu - 


9^ 

rj ■ I 


zJ 4 ^ 

A=±l 


(l)et*il)-rr,^ 


iv • ir 


( 12 ) 


is the numerator of the light-cone gauge gluon propagator with rj^ = (0,1; Oj.) in the (-1-,—,±) notation, such that 
A~^ = ry • Z = 0 is the gauge condition. For an on-shell (P = 0) gluon, the numerator of the gluon propagator in 
Eq. (1121) can be written as a sum over physical (transverse) gluon polarizations with the polarization 4-vector 


e^{l) = 



(13) 


in the light-cone gauge and = —(l/-\/2)(A,i). (Note that the last term on the right hand side of m vanishes if 
P = 0.) In arriving at Eq. (HB we have also used the fact that the quark-gluon vertex changes sign when carried 
across the cut due to complex conjugation. 

From Eq. (HID we conclude that contributions of top two diagrams in Fig. [3] can be found by only calculating the 
left diagram with the retarded gluon Green function: this conclusion is illustrated in the bottom diagram in Fig. [31 
with the retarded gluon Green function denoted by a bold arrow on the gluon line. Note also that the retarded Green 
function arising in Eq. (HID implies that > xf in the notation shown in Fig. [31 such that the gluon is first emitted 
by the quark line, and then it enters the triple gluon vertex leading to the production of the measured gluon. The 
arrow in the last diagram of Fig. [3] indicates this direction of time flow. 

So far we have considered an example of just two diagrams which combine to give us a retarded gluon Green function. 
However, the statement that by moving an end of a gluon line across the cut and by adding that contribution to the 
original diagram we would obtain a retarded Green function for that gluon is valid in general, at this classical level. 
To demonstrate this explicitly we need to consider many different types of diagrams: this is done in Appendixj^ Note 
also that our ability to use a retarded gluon Green function in calculating the amplitude should not depend on what 
goes on in the complex conjugate amplitude: while naively our argument in Fig. [3] seems to depend on the absence of 
final-state (post-shock wave) interactions in the complex conjugate amplitude, the argument is in fact true in general, 
as long as we are working at the classical level, that is, calculating order-a^ cross-section contribution involving two 
interacting projectile nucleons (quarks). We refer the reader to Appendix El for details. 

The diagram in the top left panel of Fig. |3|is in the (i)-class by the classification presented in Fig. [31 while the 
diagram in the top right panel of Fig. [3] is in the (iii)-class. We see that the contribution of the type-(iii) diagram is 
included in the type-(i) diagram by using a retarded gluon propagator in the latter. Again this conclusion is true in 
general: type-(hi) diagrams are included in the calculation of type-(i) and type-(ii) diagrams if we use retarded gluon 
propagators. 






FIG. 4. Sample diagrams with a gluon exchange between the projectile quarks. 
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Now let us consider diagrams where the projectile quarks exchange a gluon with each other (the gluon may or may 
not go through the shock wave). An example of several such diagrams is given in Fig. 01 where graphs labeled (1), 
(2) and (3) are different from each other only by the placement of the gluon exchanged between the projectile quarks. 
Once again, the arrow on this gluon line denotes a retarded gluon Green function, with the direction of the arrow 
indicating the direction of the resulting (light-cone) time ordering. We see that each diagram in Fig. |4] represents a 
sum of two diagrams (akin to Fig.jS]): each graph of Fig. 0] implies a sum of itself (with the Feynman gluon propagator 
for the gluon with the arrow on it) and the contribution where the lower quark-gluon vertex of the same gluon is 
carried across the cut (to the region before the shock wave interaction in the complex conjugate amplitude), resulting 
in the retarded Green function for the gluon marked by an arrow. Note that the time-ordering of the retarded Green 
function for the gluon in the graph of Fig. 0] again points towards the produced gluon, similar to the case illustrated 
in Fig. 01 



FIG. 5. An example of diagrams cancellation resulting from moving a retarded gluon propagator (denoted by the arrow 
indicating the direction of time-ordering) across the final-state cut. Any possible interaction may happen at a;'*' > 0 both in 
the amplitude and in the complex conjugate amplitude, as indicated by a shaded region. 


To analyze the diagrams in Fig. 0] we need the following additional observation: moving a retarded gluon Green 
function across the cut flips the sign of the contribution. This is illustrated by a simple example in Fig. 0] (see 
Appendixfor more details), where the cancellation is valid independent of the interactions which may happen after 
the shock wave interaction (i.e., for > 0 shown by the shaded region in Fig. 0|) on both sides of the cut. 

Employing the cancellation of Fig. 01 we see that diagram (1) in Fig. 0] is canceled by the contribution to diagram 
(3) coming from the region, 


(1) + (3)4>.+ =0. (14) 

Due to the difference in color factors, such cancellation does not happen between the diagram (2) in Fig. 0] and the 
diagram (3) with : instead, if we write the contribution of diagram (2) factoring out the color factor as 

(2) = M, we get 


(2) + (3)4>,+ = r, t^] M = (2) with r, (15) 

This conclusion is illustrated in the lower right diagram of Fig. 0] the sum of diagrams (1), (2) and (3) is simply 
diagram (2) with the commutator instead of the product of the fundamental generators in the amplitude. Once 
again, while we have illustrated the point by a simple example, it is valid in all cases: contributions of diagrams 
with the single or double gluon exchange between the two valence quarks can all be found by simply calculating a 
small subset of those graphs replacing products of t“’s in them by commutators. The details of how this happens 
and which diagrams survive in the end are given in Appendix]^ At order- the surviving diagrams are B 2 , Bg and 
Bii from Fig. [7] below (indeed diagram Bg is the amplitude to the left of the cut in the graph (2) of Fig. 0]), which 
have to be calculated with the commutators instead of the products of fundamental generators. As we will see below, 
B 2 = 0 (even with the commutator), such that only produced gluon emissions by projectile quarks after the shock 
wave interaction, that is for x'^ > 0, survive. This observation completes the analogy between the Feynman diagrams 
in the shock wave formalism and the diagrams representing the classical gluon fields in the MV model. 

We would like to stress that to establish this diagrams versus classical fields analogy, and simplify our calculations 
in the process, we had to utilize the fact that the amplitude we are calculating would have to be squared to obtain the 
cross section. Hence we have absorbed some contributions from the complex conjugate amplitude into our amplitude 
by employing the retarded Green functions and, for some graphs, commutators. Therefore, strictly-speaking, below 
we will not be calculating a standard Feynman-diagram amplitude, but a somewhat modified amplitude with retarded 
“propagators” and commutators: the square of this amplitude would still give the gluon production cross section. 

Let us summarize the main conclusions of this Section. First of all, to find the first saturation correction in the 
projectile, we only need to calculate 0{g^) and 0{g^) gluon production amplitudes with the retarded gluon Green 












functions. The time-ordering of the retarded “propagators” should be such that the gluon would be emitted first, 
and then would participate in an interaction ultimately leading to the production of the tagged gluon. Secondly, 
the contributions of the diagrams with the gluon exchanges between the projectile nucleons (quarks) can be more 
economically constructed by calculating a subset of those graph with the commutators of fundamental color matrices. 


III. GLUON PRODUCTION AMPLITUDE AT ORDER-g® 

A. Graphs A, B, and C 

The diagrams contributing to the gluon production amplitude at 0{g^) that include interactions with both nucleons 
are shown in Figs. in]and[71 The diagrams from Fig. [5] are labeled Ai with i = 1,..., 7 and will be referred to as the 
A-graphs. Similarly the graphs in Fig. [3are labeled Bi with i = 1,..., 12 and will be called the B-graphs. In addition 
to all the B-graphs one has to consider the diagrams similar to those in Fig. [7] but with the two projectile quarks 
interchanged: these will be labeled Ci with i = 1 ,..., 12 and referred to as C-graphs. 








FIG. 6. The diagrams containing a triple-gluon vertex. Horizontal solid lines depict two Wilson lines representing two quarks 
coming from two different nucleons in the projectile. All gluon propagators are retarded, with the time flowing in the direction 
of the produced gluon. 


To calculate the diagrams in Figs. [S] and [7] we will follow the rules established in the previous Section. All gluon 
propagators are retarded, such that time flows towards the measured gluon in the diagram. This means that the 
retarded propagator for the gluon exchanged between the two quarks in Fig. [7] is such that the lower quark-gluon 
vertex comes earlier than the upper one. The projectile quark lines are normalizes as Wilson lines (that is, each line 
is divided by 2P+ where P+ is the large “-I-” momentum in the line). 

Our ultimate goal is to obtain an amplitude in coordinate space, in a form similar to Eq. ([5]). The reason for 
that is the relative ease with which the correlators of Wilson lines are calculated in the transverse coordinate space, 
particularly in the MV model. For the diagrams in Figs. [6] and [7] we would have to perform Fourier transforms into 
transverse coordinate space, similar to the transition from Eq. © to Eq. ([HI) above. Our standard notation will be to 
label the projectile quarks 1 and 2, such that their transverse positions are 6 i_l and b 2 ^_■ Note also that a priori the 
two quarks in the a;“'"-direction moving projectile have different x~ positions, labeled here as 67 and 67 (see diagram 
Ai in Eig.[6l). Since the projectile is also a large nucleus, the differenc e &7 — b^, while suppressed by a power of energy, 
is enhanced by a power of the projectile atomic number Ai [III . [T^ . Il7j |: in the following we will keep the difference 
67 ~ b 2 non-zero throughout the calculation while remembering that it is sub-eikonally small, and will put it to zero 
at the end of the calculation where allowed. The difference 67 —1»7 serve as a regulator in the l ong itudinal Eourier 
transform. This is in accordance to the standard regularization of singularities in the MV model |lll. TT^. 

To obtain an amplitude dependent on b^ and 62 need to perform a longitudinal Eourier transform integrating 
over the “-I-” component of the momentum exchanged between the projectile quarks in the diagrams of Eigs. |6]and[71 
Since we are working in the A'^ = 0 light-cone gauge with the numerator of the gluon propagator given in Eq. m, 
we have to specify the regularization of the light-cone poles at = 0 in order to perform this longitudinal Eourier 
transform. Each way of regulating the light-cone pole is equivalent to specifying a particular sub-gauge within the 
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‘ yi. 

Bi 

B 2 


B 3 

B 4 
















v± 


Bq 

Bin 

B,i 

Bi9 


FIG. 7. The connected diagrams without a triple-gluon vertex. The remaining twelve graphs of this type are obtained by 
swapping the Wilson lines in the diagrams shown, Ci{bi, 62 ; 1, 2) = Bi{b 2 , fei; 2,1). All gluon propagators are retarded, with the 
time flowing in the direction of the produced gluon. 


light-cone gauge. In our calculation we will use the principal value (PV) prescription, where 




1 

J+ 


and 




0, PV 




(16) 


(17) 


The PV prescription corresponds to the held of the projectile nucleus moving in the x'^ direction satisfying the 
—>■ -boo) = —A±{x~ —>■ — 00 ) sub-gauge condition (see e.g. H). 

Other possible su b-gaug e choices include requiring that the Held of the projectile nucleus obeys the A_\_{x~ —>■ 
- 1-00 ) = 0 condition This corresponds to [13 




1+ - ie 


Xju J’fj, 

1+ + if. 


(18) 


if the momentum I Hows from the index ^ to the index v along the gluon line. The corresponding polarization 4-vector 
is 


e^{l) 



1+ + ie 



(19) 


One may also employ the A± (x —>■ — 00 ) = 0 sub-gauge condition, which results in reversing the sign of ie in Eqs. m 
and (fTOll . 

The reason we are going to use the PV sub-gauge is explained in detail in Appendix [BJ It turns out that while the 
A_l{x~ —>■ -boo) = 0 and A_\_{x~ —>• — 00 ) = 0 sub-gauges are very useful for many classical gluon field calculations 
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[HE EjlHlUl; in our shock-wave target setup the PV prescription is more economical. Namely, as detailed in 
Appendix |51 using the —>■ -|-oo) = 0 or —>■ —oo) = 0 gauge choices leads to the need to include new 

diagrams, which are not included in Figs. [5] and [7] (and are not part of the C-graphs). See for instance the first two 
diagrams in Fig.[TH]of Appendix I bI These diagrams were considered before in |4l| . In the shock-wave formalism these 
diagrams are probably classified as higher-order corrections to the interactions with the shock wave: namely, if one 
takes the leading-order gluon production diagrams from Fig.[Tl and considers an order-ofs correction to the interactions 
of either the quark or the gluon with the shock wave, one would then obtain an order-g^ contribution to the gluon 
production amplitude. Indeed such corrections are not going to be enhanced by an extra A-{ (in the cross section): 
however, shock wave interaction corrections to the lowest-order graphs from Fig. [T] involving a quark from another 
projectile nucleon would give an order-g^ contribution enhanced by a power of Ai, that is, they would be comparable 
to the A, B and C graphs. In Appendix iBl we show that while in the Aj_(x“ —>■ -boo) = 0 and A^{x~ —oo) = 0 
sub-gauges such diagrams are important, most of their contributions are zero in the PV sub-gauge. The remaining 
shock wave corrections contributions which are non-zero in the PV sub-gauge cancel in the amplitude squared, and, 
hence, can also be neglected in the calculation. We will therefore proceed by calculating the Fourier transform for all 
diagrams using the PV sub-gauge of the light-cone gauge. 

We performed all diagram calculations treating the projectile quark lines both as the regular high-energy quarks 
and as the light-cone Wilson lines from Eq. ([51) (only with non-zero b~). In the PV sub-gauge both results are 
identically equal. However, in the A_\_{x~ —>■ -l-oo) = 0 and A_\_{x~ —>■ —oo) = 0 sub-gauges we found that the 
Wilson-line approximation does not give the right answer for the B and C graphs. For more details we refer the 
reader to Appendix IbI 


1. Sample Diagram Calculation 


Since presenting a detailed calculation of each diagram is rather tedious, and would make the paper difficult to 
read, we will first work out one sample diagram, and then state the answers for the rest of the graphs. It appears 
that the diagram whose calculation illustrates most of the issues relevant to computing A, B and C graphs is A 2 . 
The diagrams is shown with all the momentum, coordinate, polarization, Lorentz-index and color labels in Fig. [5] 



FIG. 8. A more detailed rendering of the diagram A 2 . 


Using Fig. [Sjand treating the projectile quarks as light-cone Wilson lines we write the contribution of this diagram 
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as 


A2 = {ig f j dx+e^^i j dx+ J dy+ J d^xi d^X2 J ^ “''i ^ 

— OO —OO 0 —CO 

oo 

X f e-i<iriO-x+) ^-iq- {0-X+) ^ ^-il- (y+-0) d{k - 1)- ^_i(k-l)-(y+-0) ^ik-y+ 

J 27 r 27 r 27 r 27 r 

— OO 

^ iq-i I -(^1 I -bi I ) d q 2 I I _g^ I ) iij^-(x 2 j_-xij_) ifej_-(i'j_-X 2 j_) 


X / 


( 2^)2 


( 2^)2 


(27r)2 


—z 


—z 


( 2^)2 

—i 


21+ qT^ — +iel+ 2 (fc+ — /+) ^2 — ^ — 1+) 21+1 — + iel+ 

Y1 ^Ai*(9i) ^Ai (0 XI «a/(92) - 0 


2 (/c+ - 1+) {k - 1)- - {kjL-l±)^ + ie {k+ - 1+) ^ 

X g [{21 - k)^ g,p -{k + l)yg,y + {2k-l), g^p] e+ * (fc) 

X {2l+)U¥^2{k+-l+)Ui:^ 


( 20 ) 


The integrals over Xi and ccj, coming from the projectile quarks Wilson lines, along with the integral over y+, are 
regulated by the exponentials like to make them convergent at their respective infinities. This regularization 

is consistent with the regularization of Feynman propagators. Note that in Eq. (EOl) we are first Fourier-transforming 
the diagram into coordinate space, then integrating over x^, a:J and y+ with the proper light-cone time ordering and 
over xix and 5'2 _l. The shock wave is at x+ = 0. Due to the nature of the Fourier transform, the “-I-” components of 
momentum are conserved in all interaction vertices, including the interactions with the shock wave. Conversely, the 
” momentum component is not conserved at any of the vertices, such that e.g. {k — l)~ ^ k~ — l~, just like in 
light-cone perturbation theory (LCPT) [^, though our internal lines are not on mass shell. The outgoing gluon brings 
in a Fourier factor of ^ , where k~ = k‘]_/{2k+) since it is on mass shell. Note that k+ > 0. The gluon interaction 
with the shock wave brings in a factor of {2k+) [72^ g+'^ for a gluon line with light-cone momentum k+ and transverse 
coordinate x_l, according to the standard rules of the eikonal approximation (2^. Similarly, eikonal interaction of a 
quark with a shock wave yields {2P+) for a quark line with light-cone momentum P+ and transverse coordinate 

b±: as we have mentioned before, the factor of {2P'^) is removed by our normalization of the external quark lines. To 
simplify the calculations, we have replaced the numerators of the gluon propagators by polarization sums in Eq. (1201) : 
this is justified since (cf. Eq. (fT^ l 


Dp,{l) 


A=±l 


12 

(^./)2 


( 21 ) 


where the contribution of the “instantaneous” last term on the right vanishes in ^ 2 - Note that we are using retarded 
Green function regularization of all gluon propagators, in agreement with the discussion above and the arguments 
presented in Appendix El 
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Integrating Eq. (1201) over x'^, and y'^ and employing Eq. (El) yields 


A 2 = -'iig f 


3 xabc 




dl+ 

27r 


-ii+ {b--b-^ 


dq^ dq2 dl d{k — 1)' 


/ dq^ dq2 dl 

2tt 27r 27r 


27r 


qi — ie q 2 — ie I + {k — 1) — k —it 


^ ^ ^ ^ ^ ^ „igi±-{Sl±-bl±)+ig2±-{S2±-b2±)+il±-{^2±-Si±)+ikj^-{zj^-S2±) 


(27r)^ (27r)^ (27r)^ (27r)^ 
1 


21+ q^ — q^^ +iel+ 2 (fc+ — 1+) q 2 — q 2 ± + * e {k+ — 1+) 21+ I — l^+ iel+ 

X - 1 -V • 9i± e^* • d2± 

2 {k+ - 1+) {k - 1)- - ikx-lx)^ +ie {k+ - 1+) 

X [21 ■ el{k) exAl) ' - 1) - {k + 1) ■ ex,{k - 1) el{k) ■ ex,{l) + {2k - 1) ■ ex,{l) el{k) ■ cx^k - 1)] ■ (22) 

Since the expression in the square brackets of Eq. (|2^ is independent of the ” components of momenta, we can 
integrate over g)”, l~ and (fc — l)~ obtaining 


00 


d^fc d?l d?qi d'^q2 


„iqi±-(Si±-bi±)+iq 2 ±-(x 2 ±-b 2 ±)+il±-(x 2 ±-xi±)+ik±-(z±-X 2 ±) 


(27r)^ (27r)^ (27r)^ (27r)^ 


{l± -akx)- e^* qi± ■ q2X - PV-^ g2_L - {lx- a kx) e^* ■ qi± - PV- qi± ■ (l± - a kj_) e^* ■ g2± 

1 — a a 


9i± 92± (l± - a k±)^ 

, (23) 


where we have also simplified the expression in the square brackets using Eq. ED again and defined a = 1+/k+. 

Integration over a in Eq. (I23L while a little lengthy, can be straightforwardly performed using residues. In the spirit 
of regulating the a-integral with — b'^, we will use 6 ^ — b^ to tell us which way to close the contour for divergent 
terms in the a-integral, and drop it in the exponent afterwards. For the convergent terms in the a-integral we can 
put 6 ^ — bi to zero from the start. To be more specific, consider the following integral (not present in Eq. 


00 

/ 


^ g-iafe+ (t)2 ) py ^ _ 1 g-*A-+ (63 ) 

27r I — a 2 


Sign (62 - 61 ) « - Sign (62 -b^). 


(24) 


First of all, without the exponential containing b^ — b^ the integral is divergent, and it is regulated by the exponential. 
In the last step in Eq. (l24l) we drop — b^ after doing the integral. Indeed we can not drop b^ — b~[ within the Sign 
function, so we leave it as is. The above is the standard procedure when calculating observables in the MV model 
[ill[l3,[l3| ■ After squaring the amplitude we would have to average the resulting cross section in the projectile nucleus 
wave function, which would include integrating over all and 6 ^. This is left for future work since, as explained 
above, we are not ready to calculate the gluon production cross section having constructed only the 0{g^) amplitude 
in this paper. 

After carrying out the a-integral Eq. (|23ll reduces to 


A2 = -2zg^r^^ / d^x,d^X2U¥^Ug^ 


d^k d^l d?qi d^q2 


Jqi±-(xi±-bi±)+iq2±-(x2±-b2±)+il±-(x2±—xi±)+ik±-(z±-X2±) 


(27r)2 (27r)2 (27r)2 (27r)2 


r*A* 


qi± ' q 2 ± 


k\ l±_ — {k± ■ l±) k± 


kl 


+ e±* • giJ. 


? 2 ± • 


dlxdlx \kxxl±\ 

k± l± ■ {kj_ - l±) -l± k±- (fcj_ - l±) 


+ex* ■ Q2± 


gi± 


k± ll - l± {k^ ■ lx) 


11 


{kx-lx)^ 

+ i Sign( 6 ^ - bi)\kx X lx\qix ■ q 2 x 


-A* 


e r " • k± 


kl 


(25) 
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where we have defined 


Finally noticing that 


/l_L X l±_ - kx ly ky lx 


r>A* 


k']_ l± — {k± ■ l±) k± 


<l2± ■ 


k± l± ■ {k± - l±) -l±k±- {k± - l±) 
qi± ■ fc_L - l± {k± ■ l±) 


= -e^* xk±k±x l± 

= q2± X (fcj_ - l±) k± X 

= qi± xl± k± X l±, 


we can further simplify Eq. (l25l) to the form given below in Eq. (I28bl) . 


(26) 


(27a) 

(27b) 

(27c) 


2 . Results 


Before performing the transverse Fourier transforms, let us first list the expressions for the contributions of all the 
diagrams in Figs. [ 6 ] and 0 in transverse momentum space. All the diagram values below are given in the notation 
where the fundamental Wilson lines are moved to the left of all the t“-matrices for the same nucleon, similar to ©: 
this way all the E’s will vanish when we square the amplitude. 

Using the PV prescription to regulate the light-cone poles, the A-diagram contributions to the amplitude are 


Ai =2g^ 


(Pk (Pi 


(27r)2 (27r)2 


^ilx-{b2±-bi±)+ikx-{zx-b2±) jjad xdbc ^c\ 

S±J \ bix Jj_\b2± ;. 


e±-k±lj_-{l±-k±) 

X ^ ^ Sign (62 - ), 


fc 2 

A2 = — 2 ig^ / Pxi Px2 I 


(28a) 


^ ^ ^ ^ ^ ^ piqi±-iSi±-bi±)+iq2±-(x2±-b2±)+il±-(S2±-Si±)+ik±-{z±-X2±) 

(27r)2 (27r)2 (27r)2 (27r)2 


<ll± <llx 


-qix ■ q2± 


ei* X kx 


r.A* 


+ el • qix 


q 2 X X {kx - lx) , - ftj. X lx 


+ ex ■ q2X ■ 


ei* ■ kx 


+ iqix- q2X —p-Sign (62 - ) 


{kx-l±)^ 


P 


Sign(fc_L X lx) 


(28b) 


A, = - 2 


•3 f j 2 f d k d q d I ^ ^ ^ ^ ^ ^ ^ ^ ^ ■) xabcjjbd Tree f xd\ f xe\ 

J J (2^)2 (2^)2 (27r)2 ® ^ ^bix ^ [^2^^ ). 

.A* r dxx {kx - lx) 


— % 


2 i/j 

f d^X J 

1 

■ ( el* X 

lldl 

.V 

. ('" 

• c 


■ l± ■ q± - eT • l± 


lx - qx— ex ■ lx ■ 


•A* r - ^J-) ■ qx 1 ^A* 


(A:j.-/j.) 


{kx - Ix)^ 

e'r ■ qx 


Sign(A:_L x /_l) 


A 4 ( 6 i , 62 ; 1 , 2 ) — Apb2pi] 2 , 1 ), 
d^k dPl 


A 5 =- 2 ig'^ 


(27r)2 (27r)2 


il±-{b 2 ±-bi±)+ik±-{z±-b 2 ±) 


1 


kx lx (^-L - ^-l)' 


jabc^bd jjee y y (y y 

’’ bi± & 2 j_ V “i-L /I V J: 


er -kxki-i -Ixki-lx- {kx - lx) x kx Sign{kx x lx) 


Aq = 0 , 

At =0. 


(28c) 

(28d) 

(28e) 

(28f) 

(28g) 
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The _B-diagrain contributions are 
d?k dpi 


Bi=- ig-^ 

B2=0, 

B3=ig^ I 




(27r)2 (27r)2 


il±-{^2±—bi±)+ik±-{^±~^2±) ^'2- 

k\ {k_L - Z_l)2 


U§t Sign( 62 -- 6 r), 

(29a) 

(29b) 


dPk dPl 


• fc I 


(27r)2 (27r)2 


Jl±-{b2±-bi±)+ik±-{z±-b2±) ^-L ' 

fci (fc± - 4)2 


(rs.,rt*). (''fe.'‘),SiEn(6J-l.r). 


- 9,-„3 f ^2„ f ^]L ^1* ' gl^ ’ g2^ 


^4 = 2ig3 d y 


(27r)2 (27r)2 (27r)2 

X (4;.4^),Sign(52--6r), 


y j-ad j jbe 

^ 44 gi 7 


Bb=- ig' 


dPk dpi 


(27r)2 (27r)2 

X (‘'fe.‘'),sig.(i.j-i.r), 


oiT ■(&2J_-blJ_ )+*fcj_-(^-L-*>2J_ ) 




kl{k^-l±r "" ' 


(vr tH‘^] 

^-L hlJ_ V “l-L / ] 


Be — 0 , 
B7=0, 
Bs=0, 

Bg = 2ig^ 
Bio = - ig' 


(29c) 


(29d) 


(29e) 

(29f) 

(29g) 

(29h) 


ggrx-(b2^-bi^)+ifc^-(-r-b2^) 


(27r)2 (27r)2 

r dPk d^l 


(fc± - l±)'^ 

:*A* 


^iij_-(b2j_-feij_)+jfcj_-(?i--&2j_) ^-L ' 


^ (27r)2 (27r)2 

X [2-Sign(6^ - &C)] 


(fcx - 1±)‘^ 


jjad jjbc (y ^d^c\ jjbe (y y\ /29i) 

f)lJ_ bi± \ Oi^ b2± V "2J_ J 2 ' ' 

jjad ^bc (y ^c^d\ ^be (y ^e\ 

bi± bij_ \ oi_L Jb2± \ ^2J_ /2 


r„__2„„3 / j2 [ d^k d'^qi d^g2 4!* ’ gi-L ' g2-L 

Bii - Mg Jdy J 

X (''E-'A- 

B,2 =2i9> fjbj illJi 5k_fc 

J J (27r)2 (27r)2 (27r)2 


fci g?j. gL 


X t'i.l EEfc (ei;..''), [> - - ED]. 


The C diagrams are obtained by swapping the two projectile Wilson lines in Fig. [3 such that 

Q(6i,62;1,2) = B,(62,&i;2,1). 


(29j) 

(29k) 

(291) 

(30) 


As explained in Appendix |Al if one includes the contributions arising from moving the gluon exchanged between 
the two projectile quarks in the B-graphs across the cut, only the diagrams B 2 , Bg and Bn need to be calculated 
with the commutators of fundamental generators on the quark-1 line. Since B 2 = 0 we only need to take Bg and Bn 
from Eqs. (I29il) and (I29kl) above and replace —)► in them. The result is 



i=l 



d k d qi d q 2 ggt, , .(g, , -.j, )+iQ 2 , -(v 1 -fc 21 )+ik 1 ■(? 1 -fci i ) ’ k± qi± ■ q 2 ± dy^ ,e\ 

(27r)2 (27r)2 (27r)2 k"^ gi_L g4 ^ '2 


X / 


cdh jjad 


rbc 


Ut Ui 

Ol_L Oi± 


jbe 


Uy" - Ui 

yj- b2± 



(31) 


where the prime over the sum indicates that it is not a simple sum of the contributions in Eqs. (I29L but a sum of all 
of those graphs and the contributions with the exchanged gluon moved across the cut. An expression similar to (1311) 
can be obtained using Eq. dSOD for the sum of C graphs (which would also include contributions arising from moving 
the gluons across the cut). 
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Adding all the A, B and C graphs by using Eqs. (l28l) . (1^ and (l30l) yields 




i=l 

.3 


i=l 

j2„ 


= — 2ig / d Xi d X2 


d^k d^l d^qi d?q2 




(27r)2 (27r)2 (27r)2 (27r)2 


1 


dlx 


-qix ■ q2X 


g2 L X (/c L - 


kl 


+ e±* • qi± 


^ j^abc 

- 25 " 


Tfbd _ 

^bi_L 


jjce _ jjce 


( 


di_l 


{kx - IxY 

d\ / TB. j.eA 

J 2 


+ el ■ <12± 


qi± X l± 

n 


Sign(fc_L X l±) 


Vr 

02_L 


d^x 


d'^k dll dlq 


(27r)2 (27r)2 (27r)2 


■(xi-b2J_)+jE-(^-L-6i-L)+*fc-L-(*-L-^-L) fO-bc I ^d\ f ^e\ 

■' V ''l-L /l V '’2-L /i 


U‘-, 


bd 


r*A* 


' Tjce 


kl 


k± l± ■q± 

lUl 




(fc_L - l±) • q± • q± 




k± {k± - l±) ■ q± 


- {vii - C) 


-A* 


US 


k± l± ■ qx 


^x {kx-lx)‘^ql ‘^Ixlx 

X* ■ dx {kx - qx) ■ lx el* ■ lx 


^x {kx - IxY <?1 

^l* ■ kx {kx - qx) ■ lx\ 


kx Ixqx qx {kx-qx)^ll 215.9! kl {kx-qx)^ll 


+ 2g^ 


2 f d^k d^qi d^q2 ,iJ^^qs^-b^^)+iq2^■{x^-b2±_)+^k^■{z^-x^) qi^-q2X el*-kx , , 

X [u:i-u:i] {Vb,y)^ 


This is our final answer for the sum of A, B and C graphs in transverse momentum space. 
Fourier transforming Eq. (I32p into transverse coordinate space is straightforward. We use 


d^q 


piq±-xj_ 


9i 


= — In-- 

27t x± a 


(33) 


along with 


d g d± 

(2^)2 ql 


% XX 

27r xl 


(34) 


and 


f dlk d^l 

J (^ (^ 


gifej_.(?j_-X2J_)+iij_'(^2J_-£lJ_) 


el* X kx 


Sign{kx X lx) 


—i 

2 TT^ 


el* • {x2x - ^ij,) 

\X2X - XlXp 


(5[(ij_ - Six) X {zx - X 2 x)]- 


(35) 


Here, again, A ^ Aqcd is the IR cutoff in each individual nucleon. 
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Performing the Fourier transform of Eq. (|32|) we obtain 


12 


12 


A+ 1 : B.+ 1 : c. 


= - / d^Xi (fx2 5[{zx - Xix) X {z± - X2_l)] 


?A* 


{X2±-Xi±) Xi±-bi± X2±-b2± 


\X 2 X-Xll_r |flJ_-6lJ_p |x 2 _L- 62 _lP 


Cx • (5'l± - &1±) Zx - Xix X2X - b2± , e± ■ {X2X - b2±) Xi± - bix Zx - X2X 


|fix - b 


1 ± 


|2 \Zx-Xix\ |x2X-&2_lP |S’2±-^2J.P | 5 'l±-&l_Lp \ZX-X2x\ 


X /' 


Tjuci _ 

^XIJ_ ^bij_ 


jjce _ jjce 


+ 1 $ - "Z 




{zx — xx) xx — b 


i± 


\ZX - a;x 


Ifx - b 


1 -L| 


XX - b2X 
|xx - b 


2 ±| 


ex • {zx - bix) Zx - XX fx - ^2± ?x ' - ^ij.) - bix xx - b2x 

\zx-bix\'^ \zx-xx\'^ \xx-b2x\^ 


\zx-bix\'^ \xx-bix\'^ \xx-b 2 xr^ 


- - C) 


r-A* 


US 


{zx — Xx) Xx — b 


1_L 


XX - k 


' 2 ± 


r*A* 


[zx - b 2 x) Zx -XX XX - bix 


\zx — xxr \xx — bix\'^ \xx—b2x\'^ \zx — b2x\'^ \zx — Xx\ 


ex • {zx - b2x) XX - bix sj. - b 


' 2 ± 


4:Tt'^ 


Zx-b2x\'^ \xx-bix\'^ \xx-b2x\ 

3 


r*' {viX"), {kj’), 


{UV - ui‘‘ ) m‘ In ■ 

' ■^_L ni I 2 n- 


&2_L 


\zx - bix\ 


Zx - & 2 ± I A 

15 I dX [< -X?] /“ 


- C/S'' (C/f - Ul‘ ) "x*' Xj- Sax) 

blJ ' hn I ' 


b2± ^ 


l^± - felAP 


\zx - & 2 ±p 


Zx - &i_l| a 


r*A* 


(zx - x±) XX - bix XX - b 2 x 


\zx—Xxr |xx — fcixp \xx—b2x\^ 


Sign (62 - 6i )• 


(36) 


This is our final answer for the A, B and C graphs. 

Previously a calculation of a similar quantity, the order-g^ classical gluon production amplitude, was carried out in 
[s^ . It would be desirable to compare the two results. Unfortunately a direct comparison appears to be impossible 
at the moment: first of all, the calculation in [s^ appears to have employed a different regularization of a-integrals 
from our use of the 65 — b^ difference. The result obtained in [52| does not depend on 65 — bi ■ In addition, the 
expression derived in [s^ l in momentum space still contains an analogue of our |fcx x /x| in the denominator (see e.g. 
our Eq. (1251) above), and is not simplified to cancel out this factor. It appears that significant simplifications of the 
result in |52l | need to be performed before it could be compared with our Eq. (I32() . Furthermore, the calculation of 
appears to be done in a gauge where there are no i? or C graphs. Therefore, it is likely that the result of [13 
does not involve the calculational tricks of moving the whole gluon across the cut in diagrams B and C detailed above 
around Fig.|4]and below in Appendix 1X1 This could be another difference between the two calculations. Clearly more 
work is needed to compare our results with those in [h^ . 


B. Graphs D and E 

There is a subset of type-(i) diagrams (in the classification of Fig. [2]) that we have not calculated yet. These are the 
diagrams where one of the nucleons in the projectile is a spectator in the amplitude, that is, it does not emit any s- 
channel gluons but may interact with the shock wave, while another nucleon is a “spectator” in the complex-conjugate 
amplitude. An example of a diagram of this type representing the amplitude squared is shown in Fig. [9l 

Clearly diagrams of the type shown to the left and right of the cut in Fig. [5] contribute at order- 5 ^ to the gluon 
production amplitude and need to be resummed. (Since both projectile valence quarks participate in the interaction 
in the amplitude squared, the contribution in Fig. IH] comes in with two powers of just like squares of A, B and 
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FIG. 9. An example of a contribution to the amplitude squared, in which one of the projectile nucleons in the diagram to the 
left of the cut does not emit any s-channel gluons, while the other projectile nucleon does not emit any gluons in the complex 
conjugate amplitude. 


C graphs.) We will label such diagrams as H-graphs: they are illustrated in Fig. [TUI Let us stress that Fig. (TUI only 
contains the H-diagrams which are not zero due to color averaging: we have to take a color trace (divided by N^) 
for each projectile quark after squaring the amplitude. Taking into account that there are no emissions from quark 2 
on the other side of the cut implies that many of the diagrams in this class give zero contribution to the amplitude 
squared: those graphs are not shown. 

There are also Fl-diagrams, which are obtained by swapping the two projectile quarks (Wilson lines) in Fig. 1101 


Ei{bi, 62 ; Ij 2 ) — Di{b2, bi; 2 , 1 ). 


(37) 








FIG. 10. The “disconnected” diagrams, where one nucleon does not emit any gluons. The remaining graphs of this type are 
obtained by interchanging the quark lines in the diagrams shown, Ei(h\, 62 ) = Di{b2, &i). 


Just like for the A, B and C graphs, in summing U-graphs one has to include the contribution where the gluons 
emitted by quark 2 are either emitter, absorbed, or both emitted and absorbed on either side of the cut. This again 
results in all gluon propagators being retarded. In addition, as explained in Appendix [^ the sum of diagrams ZI 5 
and I? 6 , along with the diagrams D'^, D'^, Dq, Dq where either part of or the whole gluon loop are in the complex 
conjugate amplitude (see Fig. USD is equal to diagram Dq with the commutator of the color matrices , t^] instead 
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of t°‘ With this in mind, the values of the D graphs are 


Di = — i / d xi d X2 


f ^ ^ ^ ^ ^ ^ piqi±-{Sl±-b2±)+iq2±-{S2±-b2±)+il±-{S2±-Si±)+ik±-{^±-^2±) 


(27r)2 (27r)2 (27r)2 (27r)2 


?A» 


9m 9m 


• 92_l 


X k± ^ q2± X {k± - l±) ^ qi± X l± 




+ er • 9m 


m r ^2 + T- Sign{kx X lx) 

(k±-/_l)2 / 


— I 


rabc Tjbd rree 

J ^X2± 


fd^x J- 

1 
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(38b) 

(38c) 

(38d) 


(38e) 


Note also that color-averaging is implied: for instance, the part of the diagram Di contribution shown in Eq. (I38all 
is the one that survives a color-trace in the space of nucleon 2 after the amplitude is squared. More specifically 
when writing (t®t ^)2 in Eqs. (l38l) we only keep the part of the expression which survives the {t^t ‘^)2 —t 5’^'^/{2Nc) 
substitution. 

The sum of all the H-graphs (including the contributions from moving the gluons across the cut, as identified by 
the prime over the sum sign) is 
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(39) 


This result, along with the sum of all the E graphs, may also be directly obtained from Eq. (I5^ by “moving” the 
quark-gluon vertices on the line of quark 1 to the line of quark 2 in A, S and C graphs with an appropriate modification 
of color factors. 
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Fourier-transforming Eq. (I39|) we obtain 
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(40) 


This is our final result for the sum of the D graphs. Using Eq. (I37l) one can easily obtain the sum of the E graphs as 
well from Eq. (|40)) . Again this result can also be obtained directly from Eq. 


C. Cross-checks 

To test our main results for the order-g^ gluon production amplitude in Eqs. (|36ll and (|40)l let us run a couple of 
cross checks. 

Eirst of all, if there is no shock wave, there should be no gluon production. A simple calculation shows that if one 
puts all U = 1 and all U = 1 then 

7 12 12 6 6 

=0, ^'i?. = 0, ^'q = 0, ^'a = 0, ^'e, = 0, (41) 

i=l i=l i=l i=l i=l 


as expected. 

Another issue is gauge invariance. To test our results for gauge invariance, the order-t;^ gluon production amplitude 
can be calculated in different sub-gauges of the light-cone gauge, other than the PV sub-gauge in which the results in 
Eqs. (1551) and (HOI) were obtained. In Appendix iBlwe show that, while the calculations of the order- 5 ^ gluon production 
amplitude in the Ax(a;~ —>■ -l-oo) = 0 or Ax(a;“ —>■ — 00 ) = 0 sub-gauges are more involved than the PV sub-gauge 
calculation, the end result is the same in all three gauges for the unprimed sum of all A, B and C graphs, and for the 
unprimed sum of D and E graphs. Gauge-invariance also holds if we use retarded gluon “propagators”. Note that 
obtaining commutators of fundamental SU(Vc) generators for the B and C graphs by moving gluons across the cut 
was shown to be a legitimate trick only in the PV sub-gauge of the light-cone gauge: it is likely that to re-derive the 
result from Eq. ((551) in the Ax(a;“ —>■ -l-oo) = 0 and Ax(a:“ —>• — 00 ) = 0 sub-gauges of the light-cone gauge, the gluon 
should also be moved across the cut in the diagrams like the first two graphs in Fig. [TSl complicating the analysis. 
Hence, strictly-speaking, we have not shown that our result in Eq. (1361) (with primed sums) is gauge-invariant: while 
it may still be gauge-invariant, at the moment we can think of it as a gauge-invariant amplitude (resulting from the 
unprimed sums) with the terms added (in the PV sub-gauge) casting it in a simpler form for obtaining the gluon 
production cross section. 


IV. OUTLOOK 

This paper is the first step in our project to calculate the first saturation correction in the projectile wave function 
to the p + A classical gluon production cross section. In this paper we calculated the order- 5 ^ gluon production 
amplitude. The main results are given in transverse coordinate space in Eqs. dMl) and (1551) . Transverse momentum 
space amplitude is given in Eqs. (1321) and (1551) . These results constitute the first saturation correction to the leading- 
order gluon production amplitude in Eq. ([5]) and in Eq. respectively. 
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To complete this project and find the order-a^ contribution to the classical gluon production cross section one needs 
to calculate the order-t;® amplitude, again involving only two projectile nucleons (see Fig. HI). Similar to the I?-graphs 
above, in calculating the order-g® amplitude we will be able to simplify the color algebra by employing the fact that 
only one of the two nucleons involved emits a gluon in the complex conjugate amplitude, as shown in the diagram (ii) 
of Fig. HJ Still, this appears to be a rather challenging calculation, and it is presently left for future work. 
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Appendix A: Moving gluons across the cut 


The purpose of this Appendix is to show that (a) the inclusion of the type-(iii) diagrams from Fig. [2] into the 
calculation can be easily accomplished by using retarded gluon Green functions instead of Feynman propagators; and 
that (b) the contributions of B and C graphs (shown in Fig. [T]) to the 0{g^) amplitude along with the contributions 
of similar diagrams to the 0{g^) amplitude can all be efficiently found by calculating a subset of those graphs with 
color-commutators instead of the usual fundamental color factors. While these statements appear unrelated, they 
seem to be hard to disentangle for the B and C types of graphs. 

Let us begin with contributions to the cross section coming from the squares of the A-graphs (defined in Fig. [SI). 
It is straightforward to see by an explicit diagram analysis that the argument depicted in Fig. [3] applies in general, 
for the amplitude and the complex conjugate amplitude contributions coming from any two A-graphs. We illustrate 
this point in Fig. Illl bv considering a slightly more involved example than that shown in Fig. |3l namely by studying 
the square of diagram Ai. The argument of Eq. m applies here as well, after the following two observations. First 
of all, the interactions of the gluon that we moved across the cut in the second and third diagrams (to the left of the 
equal sign) of Fig. |TT]with the shock wave cancel, as the gluon has the same transverse coordinate on both sides of the 
cut such that U W = 1: hence we can treat this gluon as a free gluon propagator in the first graph, and as a free cut 
propagator in the second and third graphs. Second of all, in all three diagrams to the left of the equal sign of Fig. [TT] 
the interactions of quark 1 with the shock wave cancel (for a similar reason to the gluon’s, VV^ = 1), leading to the 
same color factor on this quark line, tr t^] = d“^/2, for all three graphs. Since the longitudinal coordinate integral 
ranges of the second and third graphs in Fig. |Tl] xf > and , complement each other giving independent 

xf and X 2 integrals from —00 to 0 each, we see that the sum of the three diagrams in Fig. [TT] leads to a retarded 
gluon Green function denoted by an arrow in the last diagram in Fig. 1111 just like it was demonstrated in Eq. m- 



1 


2 




FIG. 11. Diagrams contributing to the retarded Green function for the gluon propagator in the square of the Ai graph. The 
retarded gluon Green function is labeled by the arrow in the last diagram representing the net result. 


Other contributions to the amplitude squared made out of A graphs (both in the amplitude and in the complex 
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conjugate amplitude) can be analyzed in a similar manner on a case-by-case basis, leading to the same result: all 
internal gluon lines should be represented by retarded Green functions. Note that to prove the use of retarded 
propagators for the gluon lines in the diagrams like A 2 one has to also employ the fact that the “-I-” component of 
the gluon momentum is conserved in the interaction with the shock wave. 

Now we have to analyze the contributions to the amplitude squared coming from B and C graphs, along with the 
interference terms between the A and B or C diagrams. The main tool here, in addition to Eq. (HU, would be the 
cancellation of retarded Green functions from Fig. [5] (for the gluon exchanged between the quarks either before or 
after the shock wave both in the amplitude and the complex conjugate amplitude). Indeed one can proceed with the 
diagrammatic arguments, similar to how we outlined the analysis for the square of the sum of A diagrams above. 
While we performed this diagram-by-diagram check, it is rather lengthy and we are not going to present it here due 
to a very high number of diagrams one has to consider. Instead, a somewhat more compact approach to the problem 
is based on a formal argument which we present below. 

Following [53 we denote gluon fields in the amplitude by while the gluon fields in the complex conjugate 
amplitude are labeled by These fields are not the background classical fields of the target, but rather “quantum” 
fields whose Wick contractions give us the s-channel gluon propagators used in calculating diagrams A, B, C, etc. By 
analogy to ([5|) we can define the Wilson line representing quark propagation over a finite a;“'"-interval in the amplitude 

by 


r 4 

= Pexp I 4 J dx+t°- A°-~{x+,b~,b±) >. 

[ 4 


(Al) 


The corresponding Wilson line in the complex conjugate amplitude would be l4 [x ^, , 

denotes the fact that now one uses the field 4“ in the exponent of Eq. m- 
Further we define gluon propagators as contractions 


where A in the subscript 


4“(x) Aliy) ^ Df^^ix -y) = (0|T4“(x) 4^(y)|0) = 

f d^k Dy^ik) 

J (27r)^ + ie 

(A2a) 

Alix)Aliy) = Vf^,{x-y) = (0|T4“(x) 4^(2/)|0) = 

[DF,.ix-y)T 

(A2b) 

4“(x) Aliy) ^ Dliix -y) = (0|4“(x) Aliy)\0) = - ^ 

f ^ 4'’ Dy^ik) (24 Sik^) ei-k+) 

(A2c) 

4“(x) Aliy) = Dlliy - x) = (0|4“(x) 4^(y)|0), 


(A2d) 


where is given by Eq. (1161) while T and T denote time- and anti-time-ordering respectively. (As in the 

calculation leading to Eqs. (15^ and (IMl) we will use the PV sub-gauge of the light-cone gauge.) Since only the upper 
” components of the gluon fields A“ and 4“ enter the Wilson lines such as m, we will suppress the Lorentz 
indices below. That is, we will use 


— nab- 


Df{x-y) = Df — {x-y), Vf{x-y) = Vf — {x-y), D-\x - y) = D 
Note that (cf. Eq. (ITT]) 1 
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d'^k (fc) j 

y) - = _r)Z(^x-y), 


(271-) 


— iek+ 


(A3) 


(A4a) 

(A4b) 


with ZI“g((x — y) denoting the retarded gluon Green function. Indeed we are setting up the well-known Keldysh- 
Schwinger formalism [s^, Ell > which was used in small-x physics in [^, El) El] • 

Eqs. (IA4I) is a formal expression behind the cancellations like those shown in Fig. [S] To see this we concentrate 
on the Wilson lines representing quark propagators before the shock wave, both in the amplitude and in the complex 
conjugate amplitude. (For the cancellation it does not matter which interactions take place after the scattering in the 
shock wave.) The contribution of the diagrams in Fig. [5] is proportional to 


[O’ .by [ 0 ) -00]^ 0 [ 0 , -00]^ [ 0 , -00]^ 


|0), 


(A5) 
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where the ® sign indicates that the color indices of each product V are not contracted with each other, and are 
fixed at some arbitrary values. Time-ordering for ^“-fields and anti-time-ordering for ^“-fields are denoted by 
and T _4 correspondingly. For the classical gluon field diagrams we need contractions between the gluon fields 
and connecting the two projectile quarks: that is, we need contractions of the fields at (&i_l, b^) and the fields at 
( 62 ±,&^). Indeed the diagrams in Fig. [5] are examples of such contractions. Expanding Eq. (IA5I) to the first order in 

I—I 

the fields in each El V, and analyzing only the bi 62 contractions, we get 


dx+ t° 


A°- (x +,61 (a;+, 6 i , 6 i_l) 0 J dy^ {y^,b^ ,b2x) - (?/+, &2 > ^ 2 _l) | 0 ) 

— cx) —00 

0 

= -g^ J dx+ dy+ [Df{x - y) + Vf{x - y) - D^^{x - y) - D^\y - x)] 

— 00 
0 

= -52 r 0 J dx+ dy+ [Dflt[x -y)- Df^^ix - y)] = 0 . (A 6 ) 


This is the formal proof of the cancellation shown in Eig. [5] Note that before canceling, all the correlators assembled 
into the retarded Green functions. We have also employed an abbreviated notation by suppressing (&i_L,&j”) and 
(& 2±7 in the arguments of the correlators. 

A somewhat more involved calculation demonstrates that expanding Eq. (IA5I) to the second order in the fields in 

I—I 

each V and concentrating on the 61 62 contractions again, we would also get zero. This means that adding an extra 
gluon exchanged between the projectile quarks on either side of either diagram of Fig. [5] (as long as the gluon is both 
emitted and absorbed before the scattering in the shock wave) would still generate canceling diagrams. Two-gluon 
exchange is the highest order necessary in the classical field analysis m. 






FIG. 12. The sum of the diagrams complementing the square of the graph Bq and Biq. The arrows on the gluon lines denote the 
time-flow direction for the retarded propagators, while the brackets denote the color matrix commutator. Both time-orderings 
of the quark-gluon vertices on the Wilson line representing quark 2 are implied in the last two diagrams diagrams in the 
second line. The first gluon exchanged between the projectile quarks in the first diagram of the last line has the emission time 
x'^ € [0, -|-oo] with the color matrix placed to the left of the commutator in the amplitude (with similar implication for the 
last diagram). 


Let us apply the Keldysh-Schwinger formalism described above to the analysis of the squares of B diagrams. We 
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will work out the square of the diagrams Bg + Big (as defined in Fig. Ill), along with all the other similar graphs 
obtained from the square of Bg + Big by moving gluons across the cut: contributions of other diagrams can be found 
by analogy. A sample of diagrams we need to sum up is given in the first line of Fig. 1121 Namely, we fix the produced 
gluon to be emitted at light-cone time xf in the amplitude and at Xg in the complex conjugate amplitude, and sum 
up all the graphs where the other two gluons are exchanged between the two projectile quarks, as long as the gluons 
are emitted and/or absorbed after the interaction with the shock wave in the amplitude and in the complex conjugate 
amplitude (that is, for a;+ > 0). Clearly the sample diagrams in the first line of Fig. [1^ represent only a small subset 
of all the graphs that need to be included. 

Instead of summing the diagrams we notice that due to the lack of gluon emissions and/or absorptions at a:’*' < 0, 
the interactions of the projectile quarks with the shock wave cancel. Hence the contribution of all the Bg + Big 
squared type of diagrams can be written as proportional to 


(OlT^T^Itr 


„[-|-oo,a;t] Vr [-boo,a;/"] Vr ,-[a;/',0] 


X tr 


VJ [- 1 - 00 , 0 ] Vb [- 1 - 00 , 0 ] 

b2±,b^ A ^ 'a 


| 0 ) 


(A7) 


with the expression (TAtI) containing all the gluon exchanges between the projectile quarks. 
First let us rewrite Eq. (TATI) as 


(0|T^T^ |tr 
X tr 


rJ^^^_[+co, 0 | ^,.q+=o, 0 l 




| 0 ), 


(A8) 


where the parenthesis around the V factor inside the first trace are placed there just to emphasize this term. 

Once again we are interested in contractions in Eq. (IA8I) connecting the bi and bg lines. By analogy to the analysis of 
Eq. (IA5II one can conclude that two contractions between the fields in the second trace and the fields in the parenthesis 
inside the first trace cancel. We are left with contractions between [a:J, 0]^ 14^ [^2^^]a fro™ the first trace and 

the second trace, tr V/t g]^ Wj [+oo, 0]^ , along with the contractions between V/tji/", 0]^ V/Jx]'', 0]^ from 

the first trace and the same second trace, tr V/t [-|-oo, 0]^ Vb^ [+oo, 0]^ . We may also have one contraction between 
the expression in the parenthesis of Eq. (IA8I1 and the second trace, combined with another contraction either between 
10]_4 ^bi [xg, 0]_4 from the first trace and the second trace or between [a^i", 0] 14^ [a;/", 0]^ from the first 

trace and the same second trace. This is exactly the answer illustrated in the second and third lines of Fig. [121 First 

of all, since the bg quark line brings in a factor of tr [+oo, 0 ]^ Vb^ [+oo, 0 ]^ , for each contraction of A~ (bi) coming 
from anywhere in the first trace in Eq. (IA8I) with A~ (bg) from the second trace there exists a contraction of A~ (bi) with 
A~(bg); moreover, the second contraction comes in with a negative relative sign (due to hermitean conjugation in v\ 
as compared to Va), completing the original contraction to a retarded Green function, in accordance with Eq. (I A4at . 
Similarly, for each contraction between A~ (bi) with A~(bg) there exists a contraction of A~(bi) with A~(bg), which 
comes in with the relative negative sign, again generating a retarded Green function, but now via Eq. (|A4bl) . Let 
us stress again that since we are constructing gluon production cross section corresponding to the classical field, in 
this case we only need two contractions between the bi and bg lines, corresponding to the two gluons exchanged in 
the diagrams of Fig. [TH Hence our conclusions should be understood as valid for up to two gluon exchanges, but not 
necessarily beyond. The second point we need to make is that contractions with, say, [a;/', 0]^ 14^ [2^i'7 0]a 

that the gluon in the amplitude can interact with the quark bi only at light-cone times 0 < x^ < Xi, that is, to the 
left of the emitted gluon in the diagrams of Fig. [121 Moreover, since 

A r^+ nl bo-i/: t^+ nl — TTa-b\„+ 


Vb\ [a^r, 0 ]^ e Vb, [a;+, 0 ]^ = 0 ]^ 


(A9) 


the interaction is identical to that with a gluon Wilson line connecting the points (0, 6/, 6ij_) and (a:/", 6/, 6ij_) along 
the x^ light-cone direction. In passing we note that this is similar to the modification of the quark source current 
by a gluon exchange in the classical perturbative solution of Yang-Mills equations constructed in [T^ (see Eqs. (8) 


I/Jj-boo,0] 14, [+oo,0],i 


there). To conclude we see that contractions between Vj/j®!’, 0]^ I4ja;/', 0]^ and tr 

lead to the diagram to the left of the cut in the first panel of the second row of Fig. [T^] (for one contraction) and 
to the second panel of the second row in Fig. [1^ (for two contractions). The gluon propagators become retarded 
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Green functions, and the color factors are replaced by commutators. Similar arguments for contractions between 
[sj, 0]^ Vb^ [x 2 , 0]_4 and tr [+oo, 0]^ [+oo, 0]^ give the right-hand side of the first diagram in the second 

line of Fig. [1^ along with the last diagram in that line. (Let us stress that in the last two graphs in Fig. [T^] 
we implicitly include a “crossed” contribution, where the quark-gluon vertices on quark line 2 are interchanged.) 

along with another contraction 


A contraction between Vb^[x^ ,0]^ and tr [-foo, 0]^ Vbj [-foo, 0]^ 

between Vij^ [+oo, 0]^ 14^ [-l-oo, 0]^ and tr [+oo, 0]^ H 2 [+ 00 , 0]^ give the first diagram in the last row of Fig. [T2] 

(more precisely, the part of the diagram to the left of the cut), with the last graph in this last row being a mirror 
image with respect to the final-state cut. Note that these last two diagrams in Fig.fT^are zero, since the contribution 
of quark-2 line is proportional to a color trace of a commutator. We see that to calculate the contribution to the cross 
section of the square of Bg -f Bio + . •., where the ellipsis represent the multitude of other graph in this class, we only 
need to calculate the square of Bg with the retarded gluon Green function and with the color commutator instead 
of the standard fundamental color factor, along with the interference of the 0{g^) diagrams (containing a double 
commutator or two commutators, as shown in Fig. ini) with the leading-order [0{g)) gluon production amplitude, 
shown in the last four panels of Fig. [T^l We have thus completed a demonstration of the original claim that to 
facilitate the calculation one can use retarded gluon propagators, calculate only a sub-set of B (and C) graphs with 
the commutators, along with calculating only the diagrams of the (i) and (ii) types by the classification of Fig. [2j 
A complete proof of these results involves analyzing all the possible contributions to the amplitude squared. While 
still tedious, our formalism developed above makes the proof much more straightforward by reducing the vast number 
of diagrams one needs to consider to a much smaller number of Wilson line correlators like that in Eq. (TATI) . Let us 
illustrate our technique by a couple more examples. 



FIG. 13. The sum of the diagrams complementing the interference of the graphs Bg and Bio with C 5 . The arrows on the gluon 
lines denote the time-flow direction for the retarded propagators, while the brackets denote the color matrix commutator. In 
the last diagram the gluon emission time is £ [ 0 , -too]. 


First consider the interference of the Bg + Bio + ■ ■ ■ type of graphs with C 5 (the upside-down reflection of B 5 from 
Fig. [7]). This set of graphs is depicted in Fig. [T31 where we only consider gluon exchanges between the projectile 
quarks at a::+ > 0 on the either side of the cut. Noticing that the shock wave interactions with the quark 1 cancel we 
conclude that the contribution of all the graphs of the type shown in Fig. [13] is proportional to 


(0|T^T^ tr 
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We rewrite Eq. (IA10|) as 
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(All) 


If the two contractions are between eJJ-I-oo, 0]^ 14^ [-foo, 0]^ and [-boo, 0]^ Vba [+ 00 , 0]^, then they cancel just 
like they did in the analysis of Eq. (IA5|) . We are left with two options: one may have two contractions between 
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Ni"! 0]^ Vbi [ 2 ^ 1 ": 0 ]a [+ 00 , 0]^ Vb 2 [+ 00 ,0]^, giving the answer shown in the first diagram after the equal 

sign in Fig. [T31 or one contraction between 0]^ 0]^ and [+ 00 ,0]^ [+ 00 ,0]^ and another con¬ 

traction between VJ^^J-l-oo, 0]^ VbJ-boo, 0]^ and [-boo, 0]^ Vbj [-boo, 0]^. The latter case gives the last diagram in 
Fig. [T^ Interestingly the first diagram after the equal sign in Fig. [T3| is zero, since the trace of the commutator we 
get in the color space of quark 1 is zero. Note that now the final answer in Fig. [T3] is completely absorbed into the 
interference between the 0{g^) diagrams with the 0{g) diagram: hence the diagram C 5 is absorbed into the 0{g^) 
gluon production amplitude, and is not included into our 0{g^) result presented in this work. 



FIG. 14. The sum of the diagrams complementing the interference of the graphs Bg and Bio with T 5 . The arrows on the gluon 
line denote the time-flow direction for the retarded propagators, while the brackets denote the color matrix commutator. 


Our final example demonstrating our formal technique involves the interference between the same Bq -b Biq -b ... 
type of graphs with A^. It is illustrated in Fig. [TT| Similar to the above we can see that the contribution of the 
diagrams from Fig. [TT|is proportional to 


(OlT^T^/tr 


X tr 
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(A12) 


We rewrite this as 


fabo (o|T. 4 T^ <j tr 
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biJ-.fci 


X 4 .-[4:0] ,-4,0] 


tr 


4, [4:0] ^ [4,0]^ 4 ^ [+ 00 , 0] ^ [+ 00 , 0] 


| 0 ). 

(A13) 


An important difference now is that we are looking at only one contraction, corresponding to the single gluon exchange 
between the projectile quarks in Fig. [TH Similar to the above contractions between V^V’s in parenthesis cancel. 
Single contractions involving 4 [ 2 ^ 3 ": 0]^ Vhj [i^, 0]^^ in the second trace are zero: they are proportional to a trace 
of a commutator. We are left with contractions between 4 [ 2 ^^: ^]^ ^*’[ 2 ^^: ^ 1.4 ^’^d 4 [+ 0 °: 0]_4 [+ 00 ,0]^ and 

between 4 [®i "’£]4 [^i": ®]a ^i-iid 4 [+ 0 ®: 0 ]^ [+ 00 , 0 ]^: this is exactly the answer given in the last two 

diagrams of Fig. 114] 

Note that one can repeat the above argument after moving the gluon line with color b (or c) in the first graph of 
Fig. [14] across the cut, getting the same answer with two commutator terms. Adding such contributions to the final 
result from Fig. [HI would turn the propagators of the gluon lines b and c into retarded Green functions. 

The rest of the proof would be a more or less straightforward repetition of the above examples, which we verified 
but are not going to present here. In the end the propagators of all the internal gluon lines become retarded Green 
functions. The B and C graphs contribute as a subset of these graphs with color commutators. 

This Appendix presents details of the correspondence between the “tree-level” Feynman diagrams and classical 
gluon fields. It is possible that an even more compact proof of our main results can be constructed: however, at the 
moment it is not known to the authors. 
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FIG. 15. This figure illustrates the rules for calculating diagrams Us and Dg from Fig. [TO] along with the other graphs which 
need to be included. The emissions coming from quark 1 (see e.g. Fig. O in the complex conjugate amplitude are not shown 
as they are not relevant for the calculation at hand. 


Our techniques presented above can be applied to the D and E graphs as well. However, in those cases a simple 
diagrammatic summation is possible. The D diagrams which need to be added together are shown in Fig. [Tsl (The 
emissions off of quark 1, such as those in Fig. 1^1 are not shown explicitly in Fig. [TS] since they are not needed for our 
analysis.) In the end one obtains 

0^ + 0'^ + D'l + Dq + D'q Dq = T)6(with retarded propagators and with instead of t^), (Aid) 

again in agreement with the main claims of this Appendix. The proof for the E graphs is obtained by a simple swap 
of the projectile quarks I and 2. 


Appendix B: On Sub-Gauge Invariance 

In this Appendix we demonstrate that up to order the gluon production amplitude does not depend on the 
particular sub-gauge used. Using 
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2TTiS{l^), 


we can write the gluon propagator in the A^{x 


—>• -boo) = 0 sub-gauge of the light-cone gauge as 
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(Bl) 


Here we have split up the gluon propagator in the A±{x -boo) = 0 sub-gauge into two parts: one is equivalent 
to the PV gauge propagator, the other is proportional to The gluon propagator in the Ax{x~ —>■ —oo) = 0 

sub-gauge can be written in a similar form, the only difference being the sign of the S(l'^) term. The S{1^) component 
gives rise to extra terms for a given diagram when compared with the PV gauge. Due to gauge invariance of the 
amplitude it is expected that all these extra terms from all of the diagrams cancel out. This does end up being the 
case but the cancellations are not trivial. One needs to consider more diagrams than simply those shown in the main 
text as classes A, B, C, D and E. 

There are three types of new contributions in the Aj_(x“ —>■ -boo) = 0 sub-gauge calculation of gluon production 
(as compared to the PV sub-gauge calculation): most of the A-graphs change (by AAi) in the new gauge along 
with some of the D and E graphs (by ADi and AEi respectively), some of the B and C graphs acquire non-eikonal 
contributions which are still leading-order in the Aj_(x“ —>■ -boo) = 0 sub-gauge (we will refer to those contributions 
as “pinched”) and there are new diagrams which can be identified as the shock-wave interaction corrections. All of 
these contributions will be described in detail below. These contributions are necessary when dealing with both the 
A±{x~ —>■ -boo) = 0 and A±{x~ —>■ —oo) = 0 sub-gauges of the light-cone gauge. Through the rest of this section 
we will be working in the Aj_(a;“ —> -boo) = 0 sub-gauge. As will be shown, all of these new contributions can be 
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written in terms of some sort of a “gauge rotation” acting on the ordev-g single gluon emission diagrams from Fig. [T] 
Since gauge invariance is valid separately for gluon emissions from quark 1 and quark 2, we will only perform the 
calculation explicitly for gluon emission from quark 1 . 

Pinched contributions originate from terms in the B and C diagrams which, while zero in the PV sub-gauge, are 
not zero in the A^{x~ —> -boo) = 0 sub-gauge. In order to see these extra contributions one cannot just treat the 
quark lines as Wilson lines, as was done in the rest of the paper; instead one must treat them as high energy quarks 
with large and P 2 momenta. Here we examine diagram B 3 depicted in detail in Fig. 1161 to show how this pinched 
contribution arises. 



FIG. 16. Diagram B 3 , an example of a diagram with a pinched contribution. 


For the following calculation we only consider the left side of Fig. [1^1 That is, for the quark lines we only include 
terms up to the 7 ^ from the shock wave. Here we show that this diagram can be written in terms of two components, 
the PV-sub-gauge contribution which comes from treating the quark lines as Wilson lines, and a “gauge rotation”. 
The left (before the shock wave) side of the diagram gives 
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Here = (0, ,&i±) and bi^ = ( 0,62 j^ 2 _l) denote the intersection points of the projectile quarks trajectories with 

the shock wave. We also assume that the quarks are massless. We rewrite (IB2I) more compactly as 
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where we have defined the following symbols: 
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Before continuing with the remainder of the calculation it is useful to evaluate the following in the high-energy 
limit, 
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Using similar tricks for rjpT^ and IpT^ we have 


~ [l""Up{P2)\ p—, 

« [l+Up{P 2 )\ (- 1 ). 

Let us first evaluate Eq. (IB3|) by treating the quark lines as Wilson lines. This is equivalent to taking only the 

upper “-” component of the gluon propagator (lowercase a = +, /3 = +) and then taking the high energy limit, 

i.e., taking Pi^P} large. The term proportional to 5[l^) is exactly zero for a = +,/3 = +. We end up solely with 
the PV-sub-gauge contribution of the form 




dH 

w- 


( - ) ( 21 -) 


= [j+u^{Pi)]^ [-f+Up{P2)]^{ig^{P)i{P)2) 


dH 




-1 

P + ie 


PV 


PV 


(B4) 


Now lets evaluate Eq. (IB3I1 using a complete treatment of the quark lines. This allows us to accurately calculate 
the term proportional to S{1^). To do this we first evaluate the following expressions; 
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With these results in hand we evaluate Eq. dmi) arriving at 
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The first term in the square brackets in the integrand of (IB5I) corresponds to the PV part of the gluon propagator. 
Notice how this matches the result ones gets from treating the quark lines as Wilson lines, shown in Eq. dm. This 
is important because this justifies the treatment of these diagrams in the PV gauge used in the main text and in 
Appendix lAl The second term is the gauge-dependent term proportional to Notice how it also contains a 

S{1~), which embodies the “pinching” of the singularities in the U-integral; this is why we refer to such terms as 

the “pinched” contributions. This term was hidden before when we just took the “-” component of the gluon 

propagator. We can see now that the naive eikonal approximation used in arriving at Eq. (|B4p caused us to miss this 
term. 

Now that we have separated the sub-gauge-dependent part of the diagram from the PV part we can analyze these 
separately. The PV sub-gauge contribution (the first term in the square brackets of (|B5I) 1 leads to the expression for 
the diagrams given in Eq. (I29c|) of the main text. The second term in the square brackets of (IB5I) yields, after 
integration over I, 


b+UaiPl)]^ [7+Up(P2)]2 ( 

\ |Ol_L — 02J_| Ay 

We will refer to such terms as “gauge rotations”. It will be denoted graphically by a dashed line with the end points 
at the transverse positions whose difference is inside the logarithm (for Eq. (IB6I) these would be bi± and t> 2 j_). If the 
end point is attached to a quark it would give a fundamental generator, if it is attached to a gluon it would give an 
adjoint generator. Both generators have the same color. The graphical representation of the sub-gauge-dependent 
part of B 3 , algebraically given in Eq. (IB6I) . is shown in Fig. [TT] This graphical notation is used for the rest of this 
section. 

The pinched diagram contribution arises only when the gluon is exchanged between the projectile quark either 
before any other interaction takes place or when such an exchange happens after all the interactions took place. 
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FIG. 17. The pinched contribution associated with the diagram B 3 . 


Analyzing all the B-graphs in Fig. [71 we conclude that only 4 of those have the pinched contributions when they 
are calculated in the —>■ + 00 ) = 0 (or —>■ — 00 ) = 0) sub-gauge. These are i?i, B 3 , B 5 and Biq. 

The pinched contributions coming from these diagrams are denoted here as ABi, AB3, AB5 and ABiq. All of them 
are shown in Fig. 1241 below. (Indeed there are also ACi, AC3, AC5 and ACiq: these are important when checking 
gauge-invariance of the amplitude with the gluon emitted off of quark 2.) 





FIG. 18. Two left panels: examples of shock-wave correction diagrams which need to be considered in the gauges explored in 
this Appendix. Right panel: a diagram in the Ai class which is of the same order in the coupling as the two other diagrams in 
this figure. The composite object at the bottom of the diagrams is a nucleon in the target nucleus. 


Now let us focus on the shock-wave corrections. The diagrams which we have in mind are shown in the two left 
panels of Fig. [18] (cf. [4l|). These diagrams do not belong to the A, B, C, D and E diagram types considered in 
the main text: the two-gluon interaction with the target nucleon in those graphs does not leave the nucleon in a 
color-singlet state, which is indicated by the nucleon breaking up in Fig. |T8l One concludes that these diagrams are 
higher-order corrections to the Glauber-Mueller scattering in the target, which is limited by two gluons per nucleon 
|5ll |. However, his does not allow one to simply neglect these graphs, since this correction is enhanced by a power of 
A-^ (in the cross section) due to the extra projectile nucleon involved in the interaction. In the right panel of Fig.fTSl 
we show a diagram in the Ai class, with the interaction with the target limited to a single-gluon exchange with one 
of the target nucleons. Clearly all three diagrams in Fig. [TH| are of the same order in our power counting: they are 
all order-g® and have the same projectile and target nucleons participating in the interaction. Hence the first two 
graphs in Fig. 1181 along with other similar diagrams involving more multiple rescatterings in the target, are of the 
same order as the A, B, C, D and E diagrams and have to be included in the analysis. In fact one may worry why 
we did not include them in the main-text calculation in the PV sub-gauge of the light-cone gauge: below we justify 
neglecting these diagrams in the PV sub-gauge calculation performed in the main text. 

To resum the corrections of the type shown in the left two diagrams of Fig. 1181 we have to include corrections 
like this for either one of the many nucleons involved in the shock-wave interaction. We first consider scattering of 
two projectile quarks (coming from two different nucleons) on a shock-wave target. The corrections to the scattering 
on a single target nucleon are shown in Fig. [TOl and labeled Si, S2 and S3. (Our analysis will be similarly valid for 
scattering of the gluons coming from the projectile nucleons on the target.) Note that the diagram S3 in Fig. ITOl also 
has the standard eikonal contribution where the gluon exchanged between the projectile quark lines is long-lived in 
the s-channel: those types of contributions are included in the analysis of diagrams A — E in the main text and will 
be discarded here. 

In the S diagrams shown in Fig. 1191 quarks 1 and 2 are eikonal quarks traveling in the x'^ direction. Quark 3, 
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FIG. 19. The diagrams correcting the interaction of two projectile quarks (top two straight horizontal lines) with a quark 
coming from a single nucleon in the shock wave (the thick red line). 


labeled by a thick red line, originates in a nucleon from the target and thus has a large P~ momentum and is traveling 
in the x~ direction. While the evaluation of these diagrams is explicitly done for quarks, similar results exist for the 
case where quark(s) 1 and/or 2 are/is replaced by (a) gluon(s). For this calculation we will just focus on the quark 
case. The diagrams in Fig. [15] should be understood as being inside the shock wave; they can take place at any point 
in the shock wave. Since we are summing over all possible diagrams, in a given target charge distribution there will be 
diagrams where a given target nucleon has diagrams Si, S 2 and S 3 associated with it. We have to calculate and sum 
up all three diagrams associated with each nucleon in the target. As an example we calculate diagram Si explicitly. 



FIG. 20. Diagram Si, a correction to the interaction with a single nucleon in the target shock-wave. The thick red line denotes 
a quark in the target nucleon. 


This diagram Si is shown in detail in Fig. |20l As usual, we will treat the projectile quarks as the eikonal lines along 
the light-cone, while treating the target quark as a regular quark with a large component of its momentum. 
(Note that we normalize the target quark by multiplying its spinor matrix element by 1/(2P3”).) Fourier-transforming 
into coordinate space we obtain the following contribution: 
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In arriving at Eq. (lB7l) we have used the on-shell condition for the outgoing target quark, 


1 + -q+ = 



(B8) 


to eliminate the integral. Since quarks 1 and 2 are completely eikonal we replaced their spinor matrix elements 
by g°‘ ~ ~. In the target quark spinor matrix element from Eq. (lB7l) one can only have a = -I-, _L and /3 = -I-, _L 

since the gluon propagators vanish for either a = — or /3 = — due to the gauge choice. The leading contribution in 
Pg” is given by the a,/3 =_L component. Naively one would expect such contribution to be sub-eikonal, suppressed 
by a power of . As one could see below, this is indeed the case, but only until one integrates over picking the 
pole either at 1+ = 0 or at q'^ = 0. The residues at such poles generate enhancement (due to pinching of the poles in 
the product of propagators in (IB7D 1. making the end contribution leading order. Since we are interested only in this 
pinched pole contribution, below we will evaluate all the expressions keeping the Z'*' Ri 0 approximation in mind. 

Keeping the a,f3 =T component allows us to write the product of the two gluon propagators in (IB7I) as 
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Multiplying the spinor matrix element of quark 3 from (IB7I) by q±a l ±/3 from the above expression yields 
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Here we have for simplicity assumed that P^ = (0,Pg”,0j_) and defined^j_ = —qi •yi . To evaluate the spinor matrix 
elements we assumed that the spinors are chosen in the Lepage-Brodsky basis [56l | (for the ” moving quark) and 
used the table of spinor matrix elements in (see also @) with -|- -o- — to write 

Ua'iP3 + I - q)±UcriP 3 ) = SPg" ■ q± - ia l± X q±^ 

Ua'{P 3 + I - q)^±Ua{P 3 ) = PL ’ (pL “ Zj_) - ia l± X q±. 

We have also assumed that /+ ^ l± in (IBIOI) as we are only after the contribution at the /+ = 0 (or g+ = 0) pole. In 
the last step of Eq. (IBIOI) we have used the on-shell condition (IB8I) for quark 3. 

Combining Eqs. duni), dnzi) and (1591) we have 
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Integrating over Z ,g and Z'*', keeping in mind that q'^ is set by the on shell condition (IB8|) . we arrive at the final 
result 
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Using similar techniques to calculate the other two graphs in Fig. [19] we have in total 
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where we have dropped the eikonal term in S 3 labeling the remainder S 3 . 
Summing all of these terms up we arrive at 
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Notice how the polarization {a) dependence of the target quark associated with the S diagrams goes away when one 
adds the diagrams together. As we Fourier-transform Eq. (IB13I) over transverse momenta it is convenient to split up 
the resulting equation into the following four pieces, which allow for an intuitive diagrammatic interpretation: 
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FIG. 21. Diagrammatic representation of the four contributions contained in Eqs. (IB14F The short line crossing the (thick red) 
target quark in the left two diagrams indicates that the quark is on mass-shell. The dashed line represents an instantaneous 
gluon exchange. The arrow on the dashed line entering a quark-gluon vertex denotes a commutator between the fundamental 
color generators (with the color matrix of the dashed line in the first position in the commutator). 
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The diagrams representing these four expressions are shown in Fig. 1211 Notice that these satisfy the condition 

4 

Si + S2 + S3 = 5 '. 

i=l 

Let us go over the physical significance behind the four contributions in Eqs. (IB14I) . The 5— Xg (xj — xj^) term 
means all of these interactions are instantaneous, they occur at a single x^ position. (In the end of the calculation 
all the x’*' coordinates are integrated out leading to the GM exponentiation of the result.) All of the diagrams in 
Fig. are various corrections to the eikonal scattering seen in the MV model or GM approximation. Diagrams 
and S 2 consist of the eikonal scattering of the two quarks off the same target nucleon, each of them scattering in a 
classical field of the target quark, since the quark line is put on mass shell between the rescatterings (as indicated by 
a “cut” through the line). Indeed, each term in the square brackets of Eqs. (IB14al) and (IB14b|) represents a single 
t-channel gluon exchange between the projectile quarks at bi± and 62 j_ and the target quark at . The contributions 
in (IB14a|) and (IB14b|) come in with a factor of 1/2 shown explicitly in the two left diagrams of Eig. [5TJ The only 
difference between S'l and S 2 is the color factor associated with the target quark. 

Note that the contributions S'l and S 2 also occur in the PV sub-gauge of the light-cone gauge. One can see this 
by noticing that they would remain if one drops the delta-function parts of the gluon propagators in Eq. (IB7I) (and 
in a similar calculation for 5 ' 2 ). Therefore, these contributions are sub-gauge invariant. In addition, these terms 
end up canceling out when one considers the amplitude squared (since bringing either of the t-channel gluons in 
these diagrams across the cut gives an overall minus sign jdlj), meaning that these diagrams do not effect the gluon 
production cross section. This is why they are not included in the calculation in the main text leading to Eqs. dSlD 
and (HOI) . 





Q 

0 


1 1 

r ^ 

( — r 

' 0 

\ 0 

— 

8 

1 t d 

1 . c 

3 

3 


1 

i : 8 


FIG. 22. Summing up the gauge rotations for a given target nucleon distribution. 


The other two contributions, S 3 and S'^, which arise in the A±{x~ —>■ -l-oo) = 0 (and, up to a sign, in A_l(x“ —>■ 
— 00 ) = 0) sub-gauge of the light-cone gauge but not in the PV sub-gauge, correspond to one of the projectile quark 
lines scattering in a classical field of the target while being color-rotated by the field of another projectile quark (see 
e.g. da for the color-rotation terminology). The first square brackets in each of Eqs. (|B14c() and (IB14dll contain a 
single t-channel gluon exchange corresponding to the classical field of the target quark. The second pair of square 
brackets in Eqs. (IB14cl) and (IB14dl) contain the “gauge rotation” illustrated by the dashed line in Eig. using 
the notation defined in Eq. (IB 6 I) and Fig. [T71 The arrow at the end of the dashed line indicates a commutator of 
fundamental color matrices for the dashed line “gluon” and for the “true” gluon entering the quark-gluon vertex: the 
first term in the commutator contains the color matrix for the dashed line placed to the left of the color matrix of the 
“true” gluon. 


FIG. 23. Shock wave with order-^^ corrections. 


To figure out how corrections from the diagrams S 3 and S'^ contribute to the amplitude with multiple rescatterings 
we need to add up all possible gauge rotations for multiple interactions with the target nucleons. The case of the two 
projectile quarks scattering on two target nucleons is shown in Fig. [22l On the left of the equality in Fig. [22] we have 
shown two contributions arising due to the corrections like those shown in the right two diagrams of Fig. 1211 one for 
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the first nucleons and one for the second nucleon. Analyzing the sum on the left we see that the gauge rotations inside 
the shock wave cancel out and we only end up with rotations on the outside of the charge distribution. The final result 
on the right of Fig. [23] gives an order-g^ correction to the shock waves with the dashed line without an arrow defined 
as in Eq. (IB6I) and Fig. |T71 Repeating this argument for any number of target nucleons, by including both corrections 
from the right two diagrams of Fig. j^for each nucleon, we arrive at the same conclusion: the net result of all such 
corrections is equal to a diagram with a dashed line to the right of the shock wave minus the diagram with the dashed 
line to the left of the shock wave. We see that in the sub-gauge of interest each shock wave interaction is accompanied 
by the order-t;^ corrections as shown in Fig. 1231 The dashed line contributions are given by the expression in the 
parenthesis of Eq. (IB6I) , since all the x'^ coordinates are integrated out in going from S'^ and S'^ in Eq. (IB14I) to the 
contribution to the scattering amplitude. 

Using this result it is straightforward to deduce the contributions of gauge rotations to the gluon production 
amplitude. One simply has to take an order-g gluon production amplitude for the scattering of two projectile quarks 
on the target and add all possible dashed lines to it (connecting a pair of s-channel lines) immediately to the left and 
to the right of the shock wave (with the appropriate signs). The resulting corrections to the single gluon production 
amplitude are depicted in Fig. [24] below under the category (II): they are labeled AS'i. 






(Ill) 







FIG. 24. The various gauge-dependent corrections (with the produced gluon emitted by quark 1): (I) Pinched contributions; 
(II) shock wave corrections; (III) corrections to diagrams of type A and E. Each of the labels AAi, ARi, ASi, and Ai52 denote 
both the diagram above it and the sign in front of the diagram. In the PV gauge all of these contributions are zero. Notice 
how the sum of all of these is zero, demonstrating sub-gauge invariance. 


The last contributions we have to consider are the gauge dependent parts of the A, D and E diagrams. Calculating 
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these is straightforward although time-consuming, just repeating the calculation done in the main text of the paper 
but this time including the part of the gluon propagators proportional to S{1~^) (see Eq. (IB1I) 1. It is interesting to note 
that for an “instantaneous” term such as this, it does not matter whether one uses Feynman propagators or retarded 
gluon Green functions: hence our demonstration of the sub-gauge invariance also applies if one uses retarded gluon 
Green functions, as is done in the main text. 

With this in mind all of the gauge-dependent terms of diagrams A and E (corresponding to a gluon emission from 
quark 1) are shown in Fig. [Ml and are labeled AAi and AF 2 . These terms were obtained by an explicit calculation. 
Note that AAi = Ai{A_\_{x~ —>■ -l-oo) = 0 sub-gauge) — Ai(PV sub-gauge), with the same definition for AE 2 . Also 
shown in Fig. [24]are “pinched” contributions ABi. Note that in Fig. [^the labels AA^, AS^, ASt, and AE 2 include 
both the diagram above each of them and the sign in front of the diagram. All of the gauge-dependent contributions 
to Feynman diagrams at the order shown in Fig. 1241 that is, the “pinched” contributions, shock wave corrections, 
and gauge-dependent diagram corrections cancel out, preserving gauge invariance. More precisely we have 

AA; -h AB 2 + ASi = 0; 

AA'l + A52 = 0; 

AA 4 -h A 5'4 = 0 ; 

AA 5 -|- Ai?io -l- ASj = 0; 

AE 2 + AS^ = 0 ; 

AB 5 + A53 = 0 ; 

AS 3 -h ASq = 0. (B15) 

Hence our result (the sum of all the A, B and C graphs and the sum of all the D and E graphs) is independent of 
the choice of sub-gauge in the light-cone gauge. 
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